Magnetism of quantum dot clusters: A Hubbard model study 
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Magnetic properties of two and three-dimensional clusters of quantum dots are studied with 
exact diagonalization of a generalized Hubbard model. We study the weak coupling limit, where the 
electrons interact only within a quantum dot and consider cases where the second or third harmonic 
oscillator shell is partially filled. The results show that in the case of half-filled shell the magnetism 
is determined by the antiferromagnetic Heisenberg model with spin 1/2, 1 or 3/2, depending on the 
number of electrons in the open shell. For other fillings the system in most cases favors a large total 
spin, indicating a ferromagnetic coupling between the dots. 

PACS numbers: 73.21.La,75.75.-fa,71.10.-w 



I. INTRODUCTION 

Experiments have shown that cluster structure can 
show stronger form of magnetism than the bulk structure 
of the same element. This interesting behavior has been 
reported for example for Fe, Ni and Co clusters^I^EIlE] 
The size of clusters in these experiments has varied from 
few to few hundreds of atoms. Also Ga has been shown to 
exhibit paramagnetic behavior with many cluster sizes^. 
The experiments have started an intensive theoretical 
epoch in material physics, where for example the Ising 
model® , the Hubbard model^ ^ , and different fo rms of th e 
Density Functional Theory (DFT), are appliecPMliiii2l. 

Experimental breakthroughs have also been achieved 
in constru cting artificial lattices from quantum 
dotaiSElIiil The artificial quantum dot molecules 
and normal molecules (or clusters) have many dif- 
ferences, where the most important one is that the 
geometrical structure of an artificial molecule is fixed 
and the degeneracy can not be removed by Jahn- Teller 
deformation. This gives room for internal symmetry 
breaking through spontaneous magnetism, superconduc- 
tivity or superfluidity. Recent experiments have in spired 
much theoretical work| i® | i7 | i« | ii^ l ^u i ^i | ^^ l ^a i ^4 | 25 | ^ ll271^ 

For a lattice with strongly correlated particles the 
generic many-particle model is the Hubbard model, 
which has been vastly studi ed in the case of one state per 
lattice site (for reviews see^SEH). in a one-dimensional 
lattice or in a ring the Hubbard model is exactly solv- 
able using the Bethe ansatz^^l^ The magnetism of finite 
molecules^ and quantum ringa^iEIl have also been stud- 
ied using the simple Hubbard model. 

Mean-field calculations based on the spin-density func- 
tional theory, as well as ab initio configuration interaction 
calculations, predict that the Hund's first rule determines 
the total spin of an isolated quantum dot, independent 
of the interpar ticle in teraction and confining potential of 
the lattice site^ESES] ^j-^^, magnetism of the quantum 
dot molecule then depends on the total spin of the indi- 
vidual dots, on the geometry of the molecule an d on the 
coupHng between the quantum dot^SHMIllIlTlSl]^ ^ g-j^_ 
pie tight-binding model (with the exchange splitting as a 



parameter) explains qualitatively some of the magnetic 
properties of quantum dot lattices'^^ . 

Single quantum dots with a few electrons can be solved 
numerically exactly (i.e. to a high degree of conver- 
gence with respect to the necessary restrictions in Hilbert 
space) by diagonalizing the many-body Hamiltonian (for 
a review see Ref.'^^. Methods beyond the mean-field 
approxi mation have als o been applied to quantum dot 
nioleculePISE2ll3IllIlIl. 

The purpose of this paper is to study the magnetism 
of quantum dot molecules with a generalized Hubbard 
model. We consider dot molecules with two to four dots 
in two and three-dimensional geometries and up to 12 
electrons per dot. We assume the confining potential in 
each dot to be two-dimensional or three dimensional and 
harmonic at the bottom. In this case the electrons in 
one dot fill the harmonic oscillator shells. We consider 
the three lowest shells Is, Ip, and 2sld. We use a gen- 
eralization of the simple Hubbard model to describe the 
interactions: The particles interact only within a quan- 
tum dot and only the partially filled shell is included in 
the Hilbert space. The Hubbard Hamiltonian is then 
solved with exact diagonalization. 



II. THEORETICAL MODELS 



General Hubbard model 



We assume a generalized Hubbard model Hamiltonian 



H = J+U, 



(1) 



where the first term represents hopping between neigh- 
boring sites and the second term intra-site two-body in- 
teractions. 

Hoppings preserve spin, and are equal for spin-up and 
spin-down particles. Thus, J separates into two symmet- 
ric spin parts: J = J2(t=i i "^o-- The hopping part of the 
Hamiltonian is generally of the form 



nn' jj' 



nn-jj'CnjcrCn'j'aj 



(2) 



where n and n' refer to different dots, whereas j and / 
denote the orbitals in question, in n and n', respectively. 
In the case of p and d orbitals we consider directional 
dependence of the hopping integrals J (see below). In 
the case of a square shaped four-dot molecule we assume 
that there is no diagonal hopping. 

The two-body interactions of the Hubbard model are 
treated in the spirit of the tight-binding model: The elec- 
trons only interact when they are at the same quantum 
dot. Thus, U separates in the symmetric parts represent- 
ing interactions on each site n: U = X^n^n- Within a 
site, full (spin-independent) two-body interaction is con- 
sidered, which yields 



Un = 



1 



E 

J1J2J3J4 

aa' 



'~'JlJ2J3Ji'^nji(j^nj.2<T'^nJ4,a'CnJ3a (o) 



where t/ 



3l32]3]i 



are the direct space matrix elements of 
on-site interaction, depending on the interaction itself 
and the j-orbitals in question, i.e. the eigenstates of the 
confining potential. 

We do not take advantage of the fact that the Hamil- 
tonian does not depend on spin, but diagonahze the sys- 
tem for 5*2 = for even number of electrons and for 
Sz = 1/2 for odd number of electrons and only after- 
wards determine the total spin 5* for each many-particle 
state. The total number of electrons is denoted by N . 
Lanczos method was used to diagonalize the matrix. The 
largest matrix dimensions, for four dot systems with 12 
electrons, were more than 800000 with more than 1.5 
billion nonzero elements. 



B. On-site interactions 



TABLE I: Nonzero interaction matrix elements for each en- 
ergy shells of a two-dimensional harmonic confinement calcu- 
lated with the delta function interaction {U ) and with the 
Coulomb interaction {U ). Uk determines the strength of the 
contact interaction. In the case of the Coulomb interaction 
the matrix elements are given in atomic units for harmonic 
potential with eigenfrequency wq = 1 atomic units. For p- 
electrons we show for the delta function interaction the ma- 
trix elements for two representations: For px and py and for 
the angular momentum states p_ and p+. 
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We consider mainly two-dimensional systems and two 
kinds of interactions between the electrons, the delta 
function interaction and the Coulomb interaction. In 
the case of the delta function interaction and harmonic 
confinement, the relative ratios of the matrix elements 
^jiJ2J3J4. (Eq- (pi)), for each electron shell k, can be char- 
acterized by just one constant Uk- Furthermore, for the 
p-shell the ratio of the two different matrix elements is 
independent of the radial form of the confining potential. 
In the case of the Coulomb interaction we have deter- 
mined the interaction matrix elements only for the har- 
monic confinement. We consider only the open shell in 
each dot to be in the active Hilbert space and only need 
the interaction matrix elements in each shell separately. 
These are given in Table |T] Notice that for the p-shell the 
only difference in going from the short-range delta func- 
tion interaction to the long-range Coulomb interaction is 
the reduction of one U value as compared to the others. 

For triangular and tetrahedral clusters of quantum 
dots we consider also three-dimensional dots. This we 
do only within the p-shell where the nonzero matrix ele- 
ments are the same as in the two-dimensional case, but 
now also including the p^-orbitals. 



C. Hopping parameters 

In the simple Hubbard model with one state per site 
the hopping is between two s-states and independent of 
the direction. In the case of higher angular momentum 
states we have to consider also directional dependence. 
We have done this for the p-electrons and d-electrons. 

The (unnormalized) wave functions of the two- 
dimensional harmonic confinement are 



(r|l.) 
{r\p±) 
(r|2,s) 
{r\d±) 



re-"'' e=^*"^ 



= (2ar^ - l)e- 
= r2e-"'''e±2»^ 



(or (r|p5) = ^e 



(4) 



where a depends on the oscillator parameter of the con- 
finement, ^ is either x or y, and r — y'x'^ + y'^. For 
p-electrons we use both angular momentum presentation 
and Cartesian coordinates. Naturally the results are in- 
dependent of the choise of the coordinate system. We use 
the Cartesian coodinates for studying three-dimensional 



dots with p-electrons. 
and ^ can also be z. 



In that case r — ^Jx^ + y"^ + z"^ 




FIG. 1: Schematic picture of triangular and tetrahedral clus- 
ters considered. Left: triangle with pa; and Py orbitals; center: 
triangle with p^c, Py and pz orbitals; right: tetrahedron with 
Px, Py and Pz orbitals. 

Figure [T|shows triangular and tetrahedral clusters with 
p-orbitals. At the vertices electrons can occupy states 
Px, Py (and also Pz in 3D cases). We assume that the 
hopping probability is proportional to the overlap of the 
wave functions at neighboring sites. The overlap integral 
depends on the factor t — aR^, where R is the distance 
between the sites. Each nonzero overlap is proportional 
to factor 



A^ — e 
a 



-aR^/2 _ 



32 



a 



~t/2 



(5) 



Some of the overlap integrals are zero due to the sym- 
metry. The overlap integrals of the p^-electrons for a 
triangle and tetrahedron are given in Table III] In the 
case of four dots in a square the hopping between neigh- 
boring dots is along x or y direction and in this case the 
hopping from px to py is not possible. In the case of an- 
gular momentum representation the hopping probability 
also depends on the direction due the phases of the wave 
functions. Table |III| gives the hopping parameters in x 
and y directions for the p and sd shells. Note that the 
hopping probability is of the same order between all the 
states. 

We have only one parameter t describing the hopping 
probability. Notice, however, that the ratios of the differ- 
ent hopping parameters depend on t, i.e. on the distance 
between the dots. In real systems the parameter t and 
the on-site interactions U are coupled. However, we keep 
them as independent parameters in order to get a more 
general picture of the possible magnetic structures. In re- 
alistic systems where the tight-binding model (and thus 
the Hubbard model) is expected to be valid, the param- 
eter t should be of the order of 10 or larger. It turns 
out, however, that already when t > A the results are 
qualitatively independent of t. 



TABLE IL Hopping parameters Jnn'jj' (in units of A) for j 
and j' orbitals (p^ , py and pz orbitals) at vertices n and n of a 
triangle and a tetrahedron. The symbol ni refers to triangle's 
vertices 1 and 2 (see Fig. IT]), whereas symbols n2 and na go 
through all 1, 2 and 3 (n2 ^ n^). The table for tetrahedron 
shows the additional hopping parameters from the triangle to 
the fourth vertex. In the case of two signs, the upper sign is 
for ni = 1 and the lower sign for n\ = 2. 
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TABLE III: Hopping parameters for the p shell and sd shell 
calculated for the angular momentum states. J^„ijj' gives the 
asymptotic ratios (i -^ oo) of the hopping parameters (note 
that even though t — > oo, Jnn'jj' -^ due to the prefactor 
A). The hopping direction is along x and y axis. In the case 
of two signs, the upper sign is for hopping to a;-direction and 
the lower sign to y-direction. 
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Heisenberg Hamiltonian is 



1 ^ 



constant. 



(6) 



n^n' 



D. Heisenberg model 

It is well-known that in the limit of large U /t the simple 
Hubbard model with one s-state per site approaches the 
antiferromagnetic Heisenberg model. We will study if the 
clusters with half filled p or sd shells will also give the 
same spectrum as the Heisenberg model. The effective 



where L is the number of dots and Jeff is the effective 
coupling which depends on the hopping parameters and 
the on-site interaction. The spin depends on the number 
of electrons in each dot and is determined by the Hund's 
first rule. In the half-filled case the spin is 1/2, 1, and 3/2 
for the Is, Ip and 2sld shells, respectively. Note that in 
most cases of the clusters considered here the Heisenberg 
model is exactly solvable '^^ *^. 



III. RESULTS 



A. Quantum dot Dimer 
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FIG. 2; Total spin as a function of the number of electrons in 
a quantum dot dimer. The solid dots are the results for the 
lateral and vertical dimer with delta function interaction. The 
open circles and squares are the results for Coulomb interac- 
tion for lateral and vertiacal dimer, respectively. In all cases 
the ratio of the largest U and largens J is t/max/ Jmax ~ 10. 

We first consider two weakly connected quantum dots. 
Figure [2] shows the total spin as a function of electrons 
in the quantum dot dimer. In the case of the full shells, 
TV = 4, 12, 24, the spin is naturally zero. In the case 
of half filled shells the spin is also zero. However, in that 
case each dot has the spin determined by the Hund's rule, 
but the spins are opposite as in an antiferromagnet . This 
was confirmed by looking at the conditional probability. 

Figure [2] shows results for two lateral dots where the 
dots are in the same plane and for two vertical dots where 
the dots are on top of each other. In the lateral case the 
hopping probabilities are those given in Table |III| while 
in the vertical case the symmetry allows hopping is only 
between same orbitals in both dots. In the case of the 
delta function interaction the total spin is the same for 
lateral and vertical configurations. 

In general, for an open shell quantum dot the Hund's 
rule determines the spirP^. This is the case whatever 
is the the interactions between the particles, as long as 
it is repulsive. In the case of the delta function inter- 
action the favored magnetic coupling between the two 
dots is ferromagnetic. Only the half-filled case is anti- 
ferromagnetic. In the case of the long-range Coulomb 
interaction the situation is more complicated: There is a 
delicate competition between the ferromagnetic and an- 
tiferromagnetic coupling, and the total spin can also be 
between the largest possible and zero, as seen in Fig. |2] 

The antiferromagnetic order of the half-filled shell can 
be seen also by looking at the energy spectrum of the 



lowest state. The Heisenberg Hamiltonian, Eq. (l6|, for 
two spins gives energies e{S) = S{S + 1) Jcg, where S 
is the spin of the state which gets values 0, 1, • • • S'max- 
Denoting by i^s the ratio i^s = (e(S')-e(0))/(e(l)-e(0)) 
we notice that for the Heisenberg model Ks — S'(S'+l)/2. 
These ratios can now be compared to those determined 
from energy spectrum of the Hubbard model. 



TABLE IV: Ratios of energy differences Ks of half-filled dou- 
ble dots compared to those of the Heisenberg model. The 
last column gives the ratio of the largest on-site energy to the 
largest hopping integral. 
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Table |IV| shows the calculated energy ratios for differ- 
ent two dot molecules with half-filled p or sd shells. The 
results show a good agreement with the antiferomagnetic 
Heisenberg model. The results depend on the ratio of the 
on-site energies U and the hopping parameters J as indi- 
cated in the case of lateral Coulomb system. The larger 
is the ratio C/max/'/max the better the Heisenberg model 
describes the lowest energy states. 

In the case of the Coulomb interaction the spin is zero 
for all even numbers of electrons. Each dot then has an 
integer number of electrons and its spin is determined by 
the Hund's rule. The conditional probability shows that 
the coupling between the dots is antiferromagnetic. This 
result is in agreement with that of the local spin density 
approximatioEpS. 



B. Triangle of quantum dots with with p orbitals 

Next we study the magnetic properties of trimers, 
where the quantum dots form vertices of an equilateral 
triangle as shown in Fig. [T] We considered two cases: 
Two-dimensional quantum dots, where only the Px and 
Py orbitals in each dot can be occupied and a three- 
dimensional case, where each dot is spherically symmetric 
having also the Pz orbital. In the case of the triangle we 
only used the onsite matrix elements calculated with the 
delta function interaction. 

Figures [3] and |4]show results of the 2D case. In the case 
of the half-filled orbital one would expect antiferromag- 
netic coupling. Indeed, the total spin of the ground state 
is zero for A^ = 6. In the case of the triangle the antiferro- 



A 

- A, n Cj • 


Bo. 


□ o« 


i^. n o • - 



3 
2.5 

2 
1.5 

1 
0.5 





p-shell 




® 


■ 




* 


® 




@ 


® H 




n a 


- ffl D S 




a a - 


ffl D 


a 


a 



1 2 

total spin 



2 4 6 



10 12 



FIG. 3: Lowest energy levels of the Hubbard model for a tri- 
angle (without p^ electrons) with half-filled p-shell compared 
with the results of the antiferromagnetic Heisenberg model 
with spin 5=1. The distance between the dots 7? is 3 (trian- 
gles), 3.5 (squares), 4 (circles) and 5 (dots). In all cases a = 1 
and U2 = 1- The results for R = 5 agree exactly with those of 
the Heisenberg model. The energies have been scaled so that 
the ground state energy is -6. 



magnet is frustrated. The conditional probability is then 
not very helpful for confirming the antiferromagnetism. 
However, in the case of a triangle the antiferromagnetic 
Heisenberg Hamiltonian can be solved exactly and the 
energy specrum can be compared to that of the Hub- 
bard model. This is done in Fig. [3] for different values of 
the interdot distance, i.e hopping probability. The figure 
shows that all the results are in agreement with those 
of the Heisenberg model, which becomes more accurate 
when the interdot distance increases. In the case of a 
triangle the Heisenberg model has only one energy state 
for each total spin with degeneracies 1, 2, 3, and 1 for 
5 = 0, 1, 2, and 3, respectively. 

Figure |4] shows the ground state spin as a function of 
the total number of electrons in the triangle for three 
different values of the hopping parameter. In the weak 
coupling limit (black dots) the spin is at maximum for 
iV > 6, meaning ferromagnetic coupling between the 
dots. For less than half-filled case, A^ < 6, the situation is 
more complicated. The weak coupling limit (black dots) 
show maximum spin, except for N ^ 5 where S — 3/2 
(instead of the largest possible 5/2). When the couphng 
between the dots gets stronger, the total spin gets smaller 
and the tendency to ferromagnetism disappears. This is 
illustrated also in the energy plot in Fig. [4] where the en- 
ergy difference between the S = (or S — 1/2) and the 
maximum spin state is plotted for different values of the 
coupling. Weak coupling favors ferromagnetism except 
for the half-filled case while strong coupling favors small 
total spin. 

We also computed the triangle of quantum dots in- 
cluding also the Pz orbitals orthogonal to the plane of 
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FIG. 4: Upper panel: Total spin as a function of the total 
number of p-electrons in a triangular cluster of quantum dots. 
The interdot distance _R is 3 (squares), 3.5 (circles) and 4 
(dots). U2 ~ 1. Lower panel: Energy difference between 
the 5* = (or S = 1/2) state and the polarized state with 
S = N/2 calculated for C/2 = 1 and for different values of the 
hopping parameter 7J, where J is determined with R — i 
and 7 has values 0.01, 1, 4, 8, 16, 100 (from top to down). 



the triangle. Naturally, in this case each dot can occupy 
6 electrons and a half-filled dot will have spin 3/2 due to 
Hund's rule. The triangle with half-filling {N = 9) had 
again low-energy spectrum in perfect agreement with the 
Heisenberg model for spin 3/2. For other fillings the fer- 
romagnetic order was favored in the weak coupling case, 
just like in the case of the triangle without the pz or- 
bitals. It is interesting to note that the antiferromagnetic 
Heisenberg model for half-filled case seems to be indepen- 
dent of how the coupling between the individual dots is 
made. For example, when the Pz orbitals are included 
in the weak coupling limit, the hopping between them 
is nearly prohibited and their role is just to increase the 
spin per dot, i.e. the spin of the Heisenberg model. 



Tetrahedron with p^ and py and pz orbitals on 
its vertices 



The only three-dimensional cluster we considered is a 
tetrahedron. In this case we only studied the p shell 
and included all the three p-orbitals. We considered only 
the delta function interaction and used the Cartesian p- 
orbitals, the overlaps of which are given in Table In] 
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FIG. 5: Total spin and ground state energy of four quantum 
dots forming a tetrahedron, as a function of the total number 
of p-electrons (including pz electrons). The distance between 
the dots is 7? = 4 and the C/2 = 1- Note that the energy for 
the antiferromagnetic ground state for A*' = 12 is not zero but 
E = -.007533. 



The results show that for strong hopping (small dis- 
tance between the dots) and weak interaction (small U2) 
the total spin is mainly or 1 indicating no magnetism. 
For large U2 or weak coupling between the dots the re- 
sults are closely related to those observed for a triangle. 
Figure [5] shows the total ground state energy and the to- 
tal spin as a function of the number of electrons. The 
parameters correspond to the weak coupling limit. Like 
in the case of a triangle, the system is ferromagnetic for 
more than half- filling of the p shell {N > 12), while for 
less than half filling the spin is slighly reduced for some 
electron number. In exactly half-filled case the result is 
again a Heisenberg antiferromagnet with a very good ac- 
curacy. 

The total energy plotted in Fig. |5] shows that it has 
minima at iV = 5 and A'^ = 18 corresponding roughly 
1/4 and 3/4 fillings. Note that for A^ > 12 the energy 
shown is the ground state energy from which the average 
on-site repulsion, U2{N — 12), has been substracted. In 
the half- filled case the energy is only slightly negative, 
arising from 'virtual' hopping between dots. 

The asymmetry of the total spin (and energy) with re- 
spect of the half filling in the case of triangle and tetra- 
hedron can be traced back to the single particle levels of 
the tight binding model in these systems: They are not 



TABLE V: Total spin of the ground state for a four dot square 
and row and row with periodic boundary conditions (noted a 
Row-h). A' is the total number of p electrons in the four dot 
system. C refers to Coulomb interaction and S to delta func- 
tion interaction. (/max/Jmax = 17.2 for Coulomb interaction 
and 20 for delta function interaction. In the cases of square 
with 2 and 14 electrons the ground state is degenerate with 
spins and 1. 
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symmetrically distributed around the zero energy. 



D. Four dots in a square and in a row: The p-shell 

In the case of a square the hopping parameters are eas- 
ier to determine. In the case of p orbitals we have only 
hopping from px to Px and from py to Py as shown in 
the first two rows of the Table [II| for triangles. The hop- 
ping is more favorable when the orbitals point towards 
the neighboring dot. In the square this leads to an in- 
teresting situation where an electron, say in p^ orbital, 
can easily hop on the a;-direction, but can not continue 
easily around the square since every second hop would 
be in the j/-direction. The square is thus different from 
a row of four dots with periodic boundary conditions, as 
illustrated in Fig. [6] In the latter case the electron can 
continue through the row with easy hops. 




FIG. 6: Schematic picture illustrating the easy hops between 
pelectrons in a square and in a row with periodic boundary 
condition. 

Table |V] shows the total spin of the ground state of 
the four dot square as a function of the filling of the 
p shell. The result is compared to that of the row of 
four dots and the row with periodic boundary conditions. 
In the case of the Coulomb interaction the total spin of 
the of the ground state of the square of quantum dots is 
zero for all (even) numbers of p-electrons. In the case of 
the delta function interaction the cases with iV = 4 and 
iV = 12 favor maximizing the spin. The row of four dots 



favors ferromagnetism. It is interesting to note that the 
difference in the hopping probabilities between px and Py 
orbitals seem to be important for the magnetic coupling 
between the dots. For example, for iV = 6 and TV = 10 
the total spin is zero in the case of the square, while it 
has its maximum value in the case of a row with periodic 
boundary conditions. 
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FIG. 7: Energy spectrum of a square of quantum dots 
with half-filled p-levels (N=8). Black dots: results of the 
antiferromagnetic Heisenberg model. Open circles: results 
of the Hubbard model with delta function interaction and 
t/max/Jmax = 40. Plus-signs: results of the Hubbard model 
with Coulomb interaction and (7max/Jmax = 17.2. The en- 
ergies of the Heisenberg model and Coulomb interaction are 
scaled and shifted so that the lowest and highest energies 
shown agree with those of the Hubbard model with the delta 
function interaction. 

In the half-filled case each system has zero total spin. 
Figure [7]shows the comparison of the lowest energy states 
of the half-filled case with those of the Heisenberg anti- 
ferromagnetic square with S = \. For the delta func- 
tion interaction the agreement is perfect. In the case of 
the long-range Coulomb interaction the quantitave agree- 
ment is not exact but qualitatively the spectra are still 
the same. It is important to note that the directional de- 
pendence of the hopping, i.e. the fact that the electron 
can not go easily around the square without exchanging 
the orbital from px to py, does not seem to have any effect 
on the validity of the Heisenberg model. 



E. Four dots in a square and in a row: The sd-shell 

In this section wc consider four dot systems where in 
each dot the sd shell is partially filled. We assume the 
dots to be two-dimensional. The open shell then has 2s, 
d+ and d_ orbitals. The on-site interactions are given 
in Table |T] and the hopping integrals in Table |III| above. 
For simplicity we use only the asymptotic values of the 



hopping parameters and, consequently, we have only one 
free parameter which is the ratio C/max/>/max- 
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FIG. 8: 25 lowest energy states for 12 sd electrons in 
calculated with the delta function (C/max/ Jmax = 10), 
and with the Coulomb interaction (f/max/>/max = 7, 
circles, compared with the results of the Heisenber; 
black dots. The results for the Coulomb interaction 
Heisenberg model have been scaled and shifted so 
ground state and the lowest 5* = 4 state agree with 
the delta function interaction. 
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FIG. 9: 25 lowest energy states for 12 sd electrons in a row 
of four quantum dots. Black points are results for a row of 
vertical dots calculated with the delta function interaction 
(t/max/J'max = 40), opeu circles are results for a row of lateral 
dots calculated with the Coulomb interaction (f/max/Jmax = 
7.5). The results Coulomb interaction have been scaled and 
shifted so that the ground state and the lowest S = 4 state 
agree with those of the delta function interaction. 

We first studied the ground state spin as a function 
of the total (even) number of electrons for a square of 
quantum dots using the Coulomb interaction. In the case 



N = 2 and N — 22 electrons the ground state is degen- 
erate with spins and 1. For A^ = 4 the total spin of the 
ground state is 2 and for A^ = 10 it is 1. For all other 
even numbers of electrons the spin of the ground state is 
zero. In the case of a row of four lateral dots we studied 
the ground state only up to half filling. The total spin 
for the row is zero for A^ = 2, 4, 6, 8 and 12, but for 
A^ = 10 it is 3. Kolehmainen et apS studied the square 
of four quantum dots with 16 sd-electrons using the lo- 
cal spin density approximation. They found the ground 
state to be ferromagnetic which is in disagreement with 
the present result S — 0, suggesting that the LSDA can 
not describe correctly the strongly correlated system in 
this case. 

In the case of half-filling, i.e. with A^ = 12, we com- 
pared the low energy spectra with that of the Heisenberg 
model with S — 3/2. Figure IS] shows the 25 lowest en- 
ergy states for the square of four dots calculated with 
the Coulomb and delta function interactions, compared 
to the results of the Heisenberg model. The delta func- 
tion interaction reproduces nearly perfectly the spectrum 
of the Heisenberg model. Note that more levels are shown 
for the Heisenberg model and that some of the levels are 
degenerate. In the spectrum of the Coulomb interaction 
the lowest state for each spin is still in fair agreement 
with the Heisenberg model, but the higher levels do not 
agree. 

We also studied in detail the spectrum of half-filled 
row of four dots. In this case we dot not have an ana- 
lytic solution for the Heisenberg model. Instead we com- 
puted the row of four dots in the vertical arrangement us- 
ing the delta function interaction and using a large ratio 
Umax/ Jmax = 40, which we believe gives nearly exactly 
the result of the Heisenberg model (for S — 3/2). Figure 
[9]shows a comparison of that result with the spectrum of 
a lateral row of four dots with Coulomb interaction. In 
this case the agreement of with the stecra is surprisingly 
good. 



IV. CONCLUSIONS 



We have studied the magnetism of the ground sates 
of clusters of quantum dots. We described each quan- 
tum dot with a harmonic confinement and used a gener- 
alized Hubbard model for describing the interaction be- 
tween the dots. For the on-site interaction we used either 
the Coulomb interaction or the delta function interac- 
tion. Only the partially filled harmonic oscillator shell 
was used as an active basis of the many-particle calcula- 
tion. 

The quantum dot dimer favors ferromagnetic ground 
state when the interaction is short range delta function. 
Only the half-filled shells have antiferromagnetic order 
and energy spectrum in good agreement with the Heisen- 
berg model. The case of Coulomb interaction is more 
complicated, but also in this case the half-filled shells 
agree with the Heisenberg model. 

An equilateral triangle of three dots and a tetrahedron 
of four dots was studied in the case of partially filled p- 
shell. Also in these frustrated cases the half-filled cases 
are well described with the antiferromagnetic Heisenberg 
model while otherwise ferromagnetic coupling between 
the dots is favored. 

In the case of lateral four dot molecules we studied 
the row and the square. The half-filled case is again 
described well with the Heisenberg model. When the in- 
teraction is the delta function repulsion the agreement 
is nearly perfect and for Coulomb interaction also sur- 
prisingly good. This is true for both the p-shell and the 
sd-shell. For other fillings the Coulomb interaction seems 
to favor small or zero spin for the ground state. 
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